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ABSTRACT 



Devices and analytical techniques are disclosed for 
measuring the spatial profile of permittivity of a mate-, 
rial by multiple wavenumber interrogations. Electrode 
structures are disclosed which define a number of dif- 
ferent fundamental wavelengths (or wavenumbers). 
Spatially periodic interrogation signals (of temporary 
frequency "to") from the electrode structures are atten- 
uated to varying degrees by the material undergoing 
analysis, depending upon the wavenumber ("k"), 
thereby permitting the derivation of a composite dielec- 
tric profile. 

19 Claims, 4 Drawing Sheets 
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MEASURE GAIN AT SHORTEST 
YIAVELENGTH AND OERIVE 
THEREFROM A COMPLEX PERMITTIVITY 
FOR A FIRST MATERIAL LAYER 
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FIRST LAYER, DERIVE A COMPLEX PER- 
MITTIVITY FOR A SECOND MATERIAL 
LAYER 



REPEAT PROCESS 
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REFINE THE ESTIMATES OF THE FIRST AND 
SUBSEQUENT LAYERS BASED ON FIRST 
ESTIMATIONS OF ALL THE LAYERS OTHER 
THAN THE ONE BEINC ESTIMATED 
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controller. The electrode structure can be at the mate- 

APPARATUS AND M ETHODS FOR MEASURING rial surface or embedded. The electrodes are also used 

PERMITTIVITY IN MATERIALS to measure the effect of the material on a charge in- 
duced on the electrodes in response to this field. By 

BACKGROUND OF THE INVENTION 5 varying the wavenumber k, a spatial distribution of 

The technical field of this invention is dielectrometry complex permittivity is deduced as a function of the 

and, in particular, the spatial interrogation of materials temporal frequency cd, when the property depends only 

to deduce their physical properties from complex per- °n the coordinate distance, x, perpendicular to the elec- 

mittivity measurements. trode structure. 

Dielectric measurement of materials during manufac- *° Sensors, according to present invention, are useful 
turing and the like is becoming increasingly common- not only in the monitoring of polymer curing and the 
place. For example, on-line monitoring devices are now manufacture of parts which do not permit sensors to be 
available for measuring the progress of curing in parts embedded, but also in the monitoring of various other 
molded from polymeric composites. Such monitors material changes, such as the outgassing of solvents 
typically are embedded into a portion of the part which 15 from paints, the removal of moisture from coatings, the 
will be trimmed off during subsequent finishing opera- diffusion of dopants into semiconductors, and deposi- 
tions. By measuring changes in the complex dielectric tion of materials, generally. 

permittivity of the composite, the state of curing can be in one illustrated embodiment, the electrode struc- 
deduced. Such sensors can be used to control or modify ture can be formed by an array of concentrically- 
various manufacturing parameters, such as the heating 20 wound, spiraling electrodes, and the wavenumber con- 
and/or cooling rates, or the amount of applied pressure. troiier can be implemented by electronic switches 
It is known in the art that A.C. measurements of wnicn app i y voltages to subsets of the electrodes to 
dielectric properties by sensors implanted within a cur- de fine fundamental wavenumbers. In another illustrated 
ing polymer can provide useful data on curing and embodiment, the electrode array comprises a plurality 
other i Distend properties. In particular, U S.Pat No. 25 of ^ erent . gizedt terete electrode pairs (e.g., interdig- 
4,423,371 ^ued to Sentuna et al in December, 1983, itated ^ which define a lufalit of 
discloses that A.C measurements 111 1 the frequency wavenum5ers . ^ either permittivity analyzer 

range of 1 Hz to 10 kHz can be reliable indicators of - U ; ju *v • 1 * j 

curmg. See also U.S. Pat. No. 4,399,100 issued to Zsol- the resulting response sensed by the electrode 

na^aLin August, 1983, and U.S. Pat No. 4,496,697 30 "V ' T^T^ t0 P rcd ^ns based on 

issued to Zsolnly et al. in January, 1985, for further a P^^d model in order to derive a spatial profile of 

disclosures of automatic process control systems for Permittivity m the material 

curing polymeric materials. ^ mventI ° n ™U be desenbed in connection 

Conventional sensors for measuring changes in the wrth certain illustrated embodiments; however, it 

dielectric properties of a curing polymer are typically 35 s *°™a be appreciated that various modifications, addi- 

either formed as simple parallel plate capacitors, such as tions subtractions can be made by those skilled in 

those disclosed in U.S. Pat Nos. 4,399, 100 or 4,496,697 ^ witfa o ut departing from the spirit or scope of the 

or planar interdigitated capacitors, such as those dis- invention. For exampl e,Jn ^^^^J ^ o^SS^,, 

closed in U.S. Pat No. 4,423,371. The utility of these a "square spiral " oi 1 conc cntn^^ 

sensors is limited by their inability to provide enough 40 ahu wn ^rr^ouid be" clear thattheelectrode array can 

information to resolve distributions of parameters. The take various other spiraling p.e., l round spiraiyj' alter- 



deduction of inhomogeneities is sometimes attempted native geometric shapes so long as subsets or the elec- 

from temporal frequency response information, but this trodes are capable of defining different wavenumbers 

process does not give a unique distribution and cannot ^ or purposes of interrogation of the material, 

be implemented without assuming that the frequency 45 Moreover, analogous apparatus and methods can be 

dispersion of the material in question is known. employed based on measurements of inductance and the 

There exists a need for better permittivity measuring magnetic permeablility of the .materials, films, coatings, 

devices and methods. In particular, there exists a need etc. to electrically-induced y 8pgnetic fields at varying 

for- better devices and methods for conducting non- wavelengths. In a magnetic system, magnetic fields 

destructive interrogations of materials to determine 50 induce a complex inductance response in the material 

their physical properties across a spatial profile. undergoing study which can be used to deduce the 

sttmmap v nv tup TXTV^wTTrw spatial profile by analogous inductance analysis. 

SUMMARY OF THE INVENTION techniques disclosed herein can also be used to 

Devices and analytical techniques are disclosed for deduce the mechanical properties of materials or the 
measuring the spatial profile of permittivity of a mate* 55 properties of inhomogeneous electromechanical con- 
rial by multiple wavenumber interrogations. Electrode tinua, such as diffuse double layers or biphasic media., 
structures are disclosed which define a number of dif- Furthermore, the techniques for estimating the spatial 
ferent fundamental wavelengths (or wavenumbers). profile can be based on parametric models of various 
Spatially periodic interrogation signals (of temporary types, such as linear and exponential distributions of 
frequency "a>") from the electrode structures are atten- 60 permittivity, step functions, and piecewise continuous 
uated to varying degrees by the material undergoing exponentials and piecewise linear functions. In addition, 
analysis, depending upon the wavenumber ("k")i t° e techniques disclosed herein for deducing one dimen- 
thereby permitting the derivation of a composite dielec- sional spatial profiles can be extended to two and three 
trie profile. dimensional interrogations. 

The techniques described herein are referred to as an 65 ' 

"imposed o>-k" approach to dielectrometry. A spatially BRJEF DESCRIPTION OF THE DRAWINGS 

periodic field is imposed upon the material via an elec- FIG. 1 is an overall schematic diagram of a dielectric 

trode structure under the control of a wavenumber analyzing apparatus according to the present invention. 
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FIG. 2 is a schematic top view of an electrode array where the area of the first electrode is A. Measurement 
for use in the apparatus of FIG. 1. of C(a>) at a single temporal frequency therefore gives 

FIG. 3 is a schematic illustration of wavenumber at best a spatial average of the complex permittivity, 
control circuitry and resulting waveforms for use in the With assumptions concerning the frequency depen- 
apparatus of FIG. 1. 5 dence of the local complex permittivity, it is possible to 

FIG. 4 is a schematic illustration of alternative wave- use the full temporal frequency response to distinguish 
number control circuitry and resulting waveforms for between certain attributes of the non-uniformity. How- 
use in the apparatus of FIG. 1. ever, the necessity for making assumptions of this type 
FIG. 5 is a schematic top view of an alternative elec- strongly limits what additional information can be ob- 
trode array for use in the apparatus of FIG. L to tained concerning the physical mechanisms underlying 
FIG. 6 is a flow diagram of a permittivity analyzer the complex permittivity (such as molecular or ionic 
according to the present invention. dynamics). Moreover, even with these assumptions, 
___ A TT there is no unique relationship between the actual spa- 
DETAILED DESCRIPTION tjal distribution of complex permittivity and the tcm- 
In conventional dielectrometry, the temporal fre- 15 poral frequency response. This follows from Eq. 6 
quency response is used to study electrical attributes of which makes it evident that at a given frequency there 
materials, such as can be discerned through the effect on are many different distributions giving rise to the same 
the fields of molecular and domain-level polarization complex capacitance. Contributions to the integral 
and of ionic conduction (for examples). In conventional coming from the neighborhood at x=a can be inter- 
dielectric sensors, Plane Parallel electrodes and ajvolt- 20 changed with those at x=b (where a and b are any pah- 
age source, v, are used to impose an electric field, E, on of values between x=o and x=h), and the integral will 
a sample. With the charge on the driven electrode of be the same. Hence, there can be no unique inference of 
this conventional apparatus defined as q, so that the the spatial distribution from the temporal frequency 
current, i=dq/dt, the response to a sinusoidal drive response alone- 
having the frequency ct> can be defined as 25 In the present invention, devices and methods are 

disclosed for using the spatial frequency response (the 
=$/v (1) response as a function of the dominant wavenumber k) 

to deduce the spatial distribution of complex permittiv- 
where v and q are the complex amplitudes of the charge ity (as a function of temporal frequency k) when that 
and voltage. 30 property depends only on the coordinate, x, perpendic- 

In complex notation, the electroquasistatic fields in ular to the electrode structure. This approach exploits a 
the material are described by Gauss' Law property of electroquasistatic fields. If the fields are 

made to be periodic in one direction, y (the direction of 
v **~° <2> periodicity in the plane of an electrode array), they 

. 5 a/ x . A - t - A a . . 35 generally decay in the perpendicular direction, x, at a 

where 6= € (x) and the requirement that E be irrota- rate ^ h mversd proportional to the wavelength, 
tional and, hence, related to the complex amplitude of x=2 ,rA. Thus, if a spatially periodic array of elec- 
tee Potential, <P, by trodes ^ }Jse6 tQ ^p^g a potential and sense the result- 

ing charge induced on its surface, the response to the 
40 short wavelengths will reflect the permittivity of the 
Guard electrodes, having the same potential as the material m ^ immediate vicinity of the electrodes. The 
driven electrode, are often med to make E have only an f^TT ^ T ? TP* ^ l*"?* 

x component. This is done in an attempt to make the ? v '* ° f "T^X. ^ u 
electric field imposed on the sample to be uniform. A< furthe ^ °* 38 W ? U : ^ s the wavelen 8* * ^er in- 
However, the complex permittivity, y, often actually 45 crca f ed ' ^e^permittivit> rat an mcreasmg range from the 
varies with x. As a result, the electric field varies with x, sensin^g surface will influence the measurement Dis- 
and the sample cannot be subjected to a uniform field. c £ scd m ^ lowing is apparatus and method for 
The effect of the non-uniform permittivity is evident in efficicntl y reconstructing the permittivity profile from 
the measured temporal frequency response which can <A measurements made as a function of wavenumber 

be derived by writing Eqs. (2) and (3), respectively 50 I™™***)- 

» r j In one "imposed cyb" dielectrometry apparatus 10 

shown in FIG. 1, a spatially periodic field is imposed on 
^~ = o W a material 12 by means of an array of electrodes 14 at 

the surface of the material. The electrode array can be 
. d& (5) 55 formed oh a substrate of insulating material backed by a 

. * = dx highly conducting material. The sample 12, which in 

general is non-uniform in the direction, x, perpendicular 
Integration of these from a first planar electrode where to the array is shown just above the electrodes. In the 
x=o, to the second Parallel plate electrode at x=h, and apparatus shown, a subset of electrodes is driven with 
use of the conditions that the potential be zero at x=o 60 the sinusoidally varying voltage Vd by wavenumber 
and be the applied voltage, V, at x=h, gives the capaci- controller 16 while another subset is connected to high- 
tance: impedance electronics of permittivity analyzer 18 that 

measure the potential V>. Alternately, these latter elec- 
_i (6) trodes can be terminated in an arbitrary load and a 

a ! h dx 1 65 m ««urcment provided by analyzer 18 of either a sensed 

0( M ) = -J- = ^ I / voltage or current. 

L 0 J An electrode system which can be conveniently con- 

structed on a single surface and in which the wave- 
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length is varied by terminal connections is shown in trades, as well as the properties of their substrate, into 

FIG. 2. In the example shown, twelve electrodes form account. This is especially true if, as to be described 

a spiral with the direction, y, of periodicity the radial next, the inversion procedure "marches" from short to 

directions, respectively perpendicular to the four seg- long-wave lengths. Although finite difference and finite 

ments comprising one spiral. Connected for the shortest 5 element numerical methods can be used to take into 

wavelength, the electrodes have the potential shown by account the effects of geometry, a collocation technique 

switch position kj, at the top in FIG. 3. In this twelve based on a spatial Fourier decomposition of the fields 

electrode system, there are four wavelengths that can has been shown to be well suited to the direct problem, 

be obtained by appropriate connection of the electrode Using the Fourier decomposition approach, the 

teraiinations, as illustrated in FIG. 3 by k^kj. . . k* 10 mixed boundary value problem posed where the elec- 

An alternative wavenumber control system that trode array interfaces the media is solved by represent- 

achieves a variable wavelength by switching the indi- mg ^ poten tial distribution on that part of the inter- 

vidual electrode termmations is shown in Fia 4. In the face ^ md media ^ ^ betW een elec- 

sketch of the electrode-plane potential distribution lo- Xn6eg b a p iec€Mvise continuous potential having col- 

cated above the structure m FIG 4, the connections J5 locatk)n potentials> y where the ; otentials ^ <j e ced 

have been made so that ten electrodes comprise a wave- toaether 

\a ♦ * *u tr ti iT • T j specific surface capacitance density, C„ can be pre- 

are used to prevent the higher harmonics (incurred , . tK /„ rtm , qt 

because of the coupling between the edges of adjacent 20 <hcted ' ? * ?° m J> le * elec ? c .^cement re- 

electrodes) from dating the response. U otter * "™ « * e ^planeof its mterface with the 

schemes for minimizing the effects of higher harmonics, f lcctrode a complex potential in that plane 

guard electrodes can be employed at the same potential hBvm » a purely sinusoidal spatial distribution with 
as the sensor electrodes. wavenumber k* and purely sinusoidal time variation, 

An alternative embodiment is shown in FIG. 5 where „ Wlth ***** frec l uenc y *>: ™* above the elec- 
the electrode array consists of a set of different-sized trodes *** at mterface between the electrodes 
discrete pairs of electrodes. In the illustrated embodi- ( where ^ Po* 8 *^ of ^ving * surface complex per- 
ment of FIG. 5, the electrode pairs are interdigitated nuttivity, perhaps due to absorbed slightly ionized ma- 
and again formed on a substrate of insulating material terial at he mterface between the electrodes) are also 
backed by a highly conductive material. Provided the m represented by their spatial Fourier components. These 
sample facing the electrodes has the same permittivity fields m tnen made self-consistent with those in the 
profile adjacent to each pair, measurements taken from media by a numerical evaluation of the collection poten- 
each provide the desired data. Measurements of the tials, V;, and this makes it possmle to evaluate the admit- 
response of the electrode array define the gam G, which tance m ^ circuit. These admittances represent the 
in general is a complex quantity reflecting both the 35 response with arbitrary types of termination, 
magnitude and Phase of the response. In the case shown For tne particular case where the gain is defined as 
in FIG. 1, the gain might be defined as Vf/Vd- The me complex ratio of the output voltage to driving volt- 
data to be used in determining the distribution of com- a 8 e , where Yn, Yu and Y22 are the admittance of the 
plex permittivity consists of a set of gains, G=G m , "P*" network now representing the combination of 
measured at the dominant wavenumber ky, j— 1 . . . N. ^ electrode array and media subject to measurement and 

Y/ is the impedance of the circuit used to measure the 
G meari rtd=G m (u>tyf k=ki. . . Iw (7) output voltage, 

The dominant wavenumbers are related to the wave- v 

length or periodicity length, \ Jt by ^ w + ru + y/ 

ky-ZffA/ (8) _ 

The same circuit elements can be used to predict the 
The wavelength is illustrated by FIGS. 3 and 4. In the response current, as would be appropriate if the media 
following, the data designated by the subscript, j, is wcre relatively lossy over the frequency range of inter- 
taken as being in the order of increasing wavelength 50 est 

and, hence, decreasing wavenumber. Inversion methods build on having the capability for 

Various representations of the complex permittivity solving a family of direct problems. Inversion methods 
distribution can be used. For example, the distribution mu8t deduce the discrete properties needed to complete 
can be represented by layers, each having a uniform me representation from the available discrete gain mea- 
permittivity. Alternatively, each layer can have an ex- 55 surements made over an appropriate range of imposed 
ponential distribution. Included in each representation dominant wavelengths. These are the gain measurement 
is the possibility of having complex surface permittivi- summarized by Eq. 7. In general, N gains are measured 
ties at one or more of the interfaces. Without the com- ranging over wavelengths short enough to resolve the . 
plex surface Permittivities, the representations approxi- distribution close to the array and ending at a wave- 
mate the distribution by "stair-step" and piece-wise 60 length long enough to make the field extend far enough 
continuous functions, respectively. into the material to reflect material properties to the 

The "direct problem** is the basis for inferring the desired depth. With properly designed electrodes and 
distribution from the data. Here, the object is to predict material properties that do not suffer extreme varia- 
the gain, given the distribution of permittivity. To pre- tions, profiles are sensed to a depth that is approxi- 
vent errors in deducing the distribution that are associ- 65 mately a quarter wavelength. In general, the inversion 
ated with the geometry of the electrode array, h is method then provides for an identification of N parame- 
desirable to have a method of solving the direct prob- ters describing the permittivity distribution from these 
lem that takes the finite width and spacing of the elec- N measurements. In the case of a stair-step approxima- 
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tion, the distribution of complex permittivity would be refined by passing through the data, again from short to 

represented by the N complex permittivities of the long wavelengths. In these subsequent iterations, the 

steps. In the case where (perhaps exponential) functions parameters of layers other than the one being estimated 

are pieced together in a continuous manner, the N iden- are taken as being their values on the last estimation, 

titled permittivities could be those at the locations 5 An optimal parameter estimation scheme takes ad- 

where the smooth distribution is pieced together. In any vantage of ^ a priori k^^ge 0 f the complex per- 

case, the introduction of complex surface j^ttivities ^vity distribution. For example, if it is known that 

imposes a requirement for additional data. Thus, with N the ^riil actually consists of one or more discrete 

measurements over an appropriate range of wave- j acrQSS the boundaries of which it ted that 

lengths, it would be possible to deduce N-S buUcparam- io ^ ta permittivity suffers ^ ab t than 

eters descnbmg either the stair-step or the piece- wise „ ? • u . , H ; 

distributions and S complex surface permittivities. 8 model be used m which one or more of the layers 

It is important to note that the inversion process can re P r «*°* *f boundaries between actual Physical lay- 
be carried out using measurements taken at the same ers * U ,ocatlon of such abrupt changes are not 
frequency. Thus, the permittivity distributions can be 15 known » parameter estimation scheme can then be 
deduced at each frequency. The information inherent to one m which the loca tion of such a boundary is one of 
the temporal frequency response is therefore preserved f* 10 Parameters being estimated. In some instances, the 
and can be exploited for further parameter identification inversion can best be reduced to one or perhaps a se- 
purposes. (An example would be to distinguish between quence of multiple parameter estimations. For one 
polarization mechanisms and between polarization and 20 s ^ ec * m tnc art » the single parameter estimation tech- 
ion migration phenomena as they occur at various niques desribed herein are readily generalized to the 
depths in the media.) simultaneous estimation of multiple parameters. 

The estimation of parameters from experimental data We claim: 

can be formalized by defining an error that is the differ- 1. An apparatus for measuring the dielectric Proper- 

ence between the measured gains and those predicted 25 ties of a material, the apparatus comprising: 

by the direct method: an array of electrodes capable of imposing an electric 

potential in a material and sensing a resulting elec- 

*=Gnwisind-G{V)pndicttd (io) tr i ca i response; 

xi™ a • n * nC « » a wavenumber controller connected to said electrode 

the different gains txodes and thereby defining fundamental wave- 

If, as in the illustration used here, there are as many ™*f re ™* imposed electric po- 

measured gains as parameters to be determined, the tentml m . the ™ tenal; md 

estimation process then amounts to finding those values a Permittivity analyzer connected to said wavenum- 
of 9 that make 35 ber controller, including means for comparing the 

resulting responses sensed by said electrode array 
e<0)=o (ll) at different wavenumbers to predictions in order to 

derive a spatial profile of permittivity in the mate- 
If there are more measurements than parameters so that rial 

the identification problem is overspecified, the estima- 40 2. The apparatus of claim 1 where in the array of 
tion process can be represented as a minimization of a electrodes comprises a plurality of concentrically- 
quadratic norm based on these errors and possibly a wound, spiraling electrodes. 

weighting ruction. 3. The apparatus of claim 1 wherein the array of 

As shown in FIG. 6, the estimation problem can be electrodes comprises a plurality of discrete electrode 
reduced to a sequence of single parameter estimations 45 pairs defining a plurality of fundamental wavenumbers. 

{ ? -rt? 8 tW ° pafa T er estunat i 0 ?? ) l y ^ ex - 4. The apparatus of 3 wherein the electrode pairs are 

plating the properties of quasistatic fields. For example, interdigitated electrode pairs. 

oT?™ ^°,KL" r° 0th i ^ thC & 5 « ™ e apparatus of claim 1 wherein the wavenumber 

o m (ki)measured at the shortest wavelength is used in » ,. „ . 

conjunction with a single layer model to deduce the 50 c C °* tr ° Uer f ^comprises a Plurality of electronic 

complex permittivity oflhe material nearest to the eleo Pitches for applymg at least one voltage to a subset of 

trodes. With this parameter in hand, the complex per- said electrodes of said array. 

mittivity of the material somewhat further out is deter- A 6 ' . e apparatus of claim 5 wherein the switches 

mined from the gain G m (k2> measured at the next longer rurther include means for applying two distinct voltages 

wavelength. Like the first step, this one amounts to a 55 10 ^P^ate subsets of the electrodes to define a funda- 

single parameter estimation of the complex Permittivity mentaJ wavenumber. 

of the second layer in a two-layer model for the media. 7 * apparatus of claim 6 wherein the switches 

Following this procedure, the measured gain G m (k/) is further include means for defining floating electrodes, 

used in conjunction with a model representing the com- The apparatus of claim 1 where the permittivity ■ 

plex permittivity by j parameters, the last of which is 60 analyzer is a microcomputer programmed to measure 

associated with the layer furthest from the electrodes the permittivity of the material at multiple wavelengths 

(and about a quarter wavelength from the electrodes, based upon, a parametric model of the material and 

based on the wavenumber of that measurement) to further includes means for performing iterative analyses 

make a single parameter estimation. Once the first pass to obtain a refined spatial profile, 

has been made through the data, working from short to 65 9. A method for measuring the dielectric properties 

long wavelengths and, hence, performing a sequence of of a material, the method comprising: 

single parameter estimations of the profile working disposing ah array of electrodes in contact with a 

outward from the electrodes, the parameters can be material, the electrodes being adapted to impose an 
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electric potential in the material and sense a result- 
ing electric response; 

interrogating the material at multiple wavelengths 
defined by subsets of said electrodes; and 5 

sensing the resulting electrical response at multiple 
wavelengths to derive a spatial profile of at least 
one property of the material. 

10. The method of claim 9 wherein the step of inter- 10 
rogating the material further includes sequentially inter- 
rogating the material at multiple wavelengths. 

11. The method of claim 9 wherein the step of inter- 
rogating the material further includes simultaneously J5 
interrogating the material at multiple wavelengths. 

12. The method of claim 9 wherein the step of sensing 
the electrical response to derive a spatial profile further 
includes 

20 

measuring the gain at the shortest wavelength to 
derive a profile estimation for a first material layer; 
and 



,690 

10 

measuring the gain with at least one longer wave- 
length to derive a profile estimation for a subse- 
quent material layer. 

13. The method of claim 12 wherein the method 
further includes employing a parametric model to ob- 
tain said estimations for the first and subsequent mater 
rial layers. 

14. The method of claim 13 wherein a linear model is 
employed to represent the layers. 

15. The method of claim 13 wherein an exponential 
model is employed to represent the layers. 

16. The method of claim 13 wherein a continuous 
model is employed to represent the layers. 

17. The method of claim 13 wherein a step function 
model is employed to represent the layers. 

18. The method of claim 13 wherein the method 
further includes refining the estimate of each layer 
based on a prior estimation of all the layers other than 
the one being refined. 

19. The method of claim 18 wherein the step of refin- 
ing the estimate is iteratively repeated for each layer 
until a best fit to a parametric model is obtained. 
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